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RADAR STUDIES IN THE SOLAR SYSTEM 


This project entails a collaboration involving SAO, NAIC (Arecibo and Cornell), 
and JPL, including D.B. Campbell, J.F. Chandler, J.K. Harmon, S.J. Ostro, 

I. I. Shapiro, and M . A. Slade, among others. 

Goals 

We aid in a study of the solar system by means of ground-based radar. We have 
concentrated on (i) developing the ephemerides needed to acquire radar data at 
Arecibo Observatory and (ii) analyzing the resultant data to: test fundamental 
laws of gravitation; determine the size, shape, topography, and spin vectors of 
the targets; and study the surface properties of these objects, through their 
scattering law and polarization characteristics. 

Progress and Accomplishments 

We are engaged in radar observations of solar-system objects, including 
planets, natural satellites, asteroids, and comets, both as systematically 
planned targets and as "targets of opportunity." Unfortunately, the 
opportunities at Arecibo have been nil for the past year, while the radar 
system was being replaced as part of the Arecibo upgrade project, now nearly 
complete. We have taken advantage of this time to improve the software package 
for Arecibo on-site ephemeris correction based on late-arriving optical data or 
newly-available radar observations of the target. This package gives Arecibo 
the capability to generate a corrected radar ephemeris quickly on-site if 
valuable new data become available during the course of an observing run (which 
may happen for newly-discovered Earth-approaching objects). We have also 
continued the analysis of radar data and prepared articles for publication in 
collaboration with our colleagues. Several papers were published, describing 
the results of asteroid observations. 

Anticipated Results 

We plan to continue our activities in this field, both by obtaining radar 
observations of asteroids, comets, planets, and satellites and by analyzing the 
data. As a result of the upgrade, now expected to be complete by the end of 
1996, the Arecibo radar will offer a greatly improved sensitivity to all of the 
radar targets and the possibility of observing targets not previously 
detectable . 
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Radar Images of Asteroid 4179 Toutatis 

Steven J. Ostro,* R. Scott Hudson, Raymond F. Jurgens, 
Keith D. Rosema, Donald B. Campbell, Donald K. Yeomans, 
John F. Chandler, Jon D. Giorgini, Ron Winkler, Randy Rose, 

S. Denise Howard, Martin A. Slade, Phil Perillat, Irwin I. Shapiro 

Delay-Doppler images of the Earth-crossing asteroid 41 79 Toutatis achieve resolutions 
as fine as 125 nanoseconds (19 meters in range) and 8.3 millihertz (0.15 millimeter per 
second in radial velocity) and place hundreds to thousands of pixels on the asteroid, which 
appears to be several kilometers long, topographically bifurcated, and heavily cratered. 
The image sequence reveals Toutatis to be in an extremely slow, non-principal axis 
rotation state. 


Optical investigations of Earth-orbit- 
crossing asteroids (ECAs) are limited by 
their small angular sizes (typically «1 arc 
sec), but radar measurements of the distri- 
bution of echo power in time delay and 
Doppler frequency can be used to synthesize 
images of these objects. Toutatis’s approach 
to within 0.024 astronomical unit (AU) 
(9.4 lunar distances) on 8 December 1992 
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was an unprecedented opportunity for radar 
investigation of a small body. We present^ 
images that achieve areal resolutions *—100 
times finer than previously obtained (1 ) for 
any ECA and reveal an object with startling 
characteristics. 

The images were obtained with the 
Goldstone Radar in California daily from 2 
to 18 December and at the Arecibo Obser- 
vatory in Puerto Rico daily from 14 to 19 
December (Table 1). We used binary-phase 
coded wave forms (2) that provided a time 
resolution At equal to the temporal extent 
of each code element, a range resolution 
equal to At times half the speed of light, and 
a frequency resolution 

A / = 1/(RP X N c oh x ^fft) = 1/Tcoh 

( 1 ) 

where the code repetition period RP equals 
At times the number of elements in the 
code (255 for Goldstone, 8191 for Arecibo) 
and N rOH isjhe number of successive, RP- 
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long time series of decoded voltage samples 
that are coherently summed before spectral 
analysis with a fast Fourier transform of 
length Lj.f.-y (256 tor Goldstone, 2048 for 
Areciho). A power spectral estimate tor 
each delay “bin” is obtained from coherent 
analysis of T c - ol j seconds of data; we refer to 
such an image as one ‘‘look.’ 1 Spectral 
estimation noise (3) is chi-square distrib- 
uted and the fractional root-mean-square 
(rms) fluctuation in a sum of N looks is 
Thus, fine frequency resolution re- 
quires a long coherence time and 

low noise requires incoherent summation 
of many looks. However, long integrations 
are undesirable if the resolution is degrad- 
ed by smearing from either the target’s 
rotation or its translational motion with 
respect to the delay-Doppler prediction 
ephemeris, which is used to drift the sam- 
pling time base to maintain coregistration 
of echoes from any given range bin and to 
tune the receiver to the Doppler-shifted 
echo. The Doppler-prediction error Av. ph 
defines the rate of change of the delay- 
prediction error At h and hence the rate 
of delay smearing 

dAr^dt) _ —A v, rh (t) 

' dt~ - f t ; " 

where F TX is the transmitter frequency. Our 
observations progressed through several 
generations of increasingly accurate ephe- 
merides (Table 1). 

Figure 1 shows our Goldstone 0. 5- juls 
images and Areciho 0.2-pis images (4). 
These images are very different from opti- 
cal pictures. Constant-Doppler planes are 
parallel to both the line of sight (LOS) 
and the target’s apparent spin vector, and 
constant-delay planes are normal to the 
LOS. These two orthogonal sets of planes 
cut the target into rectangular cells, and 
each image pixel shows the sum of echo 
power in one cell. Beyond the leading 
edge, each cell can capture echoes from 
surface regions on either side of the appar- 
ent equator, so rhe surface-to-image map- 
ping is potentially two-to-one or even 
many-to-one. Moreover, the kilometers- 
per-hert: conversion factor for any image 
depends on the asteroid’s apparent spin 
vector, as follows. For a rigid body whose 
center of mass (COM) is at a constant 
distance from the radar, the Doppler fre- 
quency of the echo from a point r with 
respect to the COM is 


(W lCL X rhe 

(A/2) 


(3) 


where the unit vector e points from the 
COM to the radar and we ignore parallax 
and second-order terms. The target’s ap- 
parent spin vector W )pp satisfies 


W, ipr = W + \\C V (4) 

where W is rhe intrinsic spin vector and 
W sk is the contribution from the target’s 
plane-of-sky (POS) motion (Table 1). The 
kilomctcrs-per-hert: conversion factors are 

Goldstone: 87.2/ 1 W il)T X el 

= P/(41 3cos5) (5) 

Areciho: 512/1 W, l[T X el 

= P/( l . 1 5cos5) (6) 

where spin rates are in degrees per day, P 
= 560/1 W I is the instantaneous, appar- 
ent spin period in days, and 5 — 

cos 1 ( I W „ r X e I / I W ipp I ) i> the instan- 
taneous, apparent, target-centered decli- 
nation of the radar. 

The following discussion focuses on the 
Goldstone images. The distribution ot 
echo power is extremely himodal in most 
images (two lobes were perceptible even 
on 5, 7, 14, and 17 December, when they 
were at similar ranges and overlap), indi- 
cating a prominent bifurcation in Touta- 


1 1 s’s >hnpe. On 5 to 6 and 15 to 16 De- 
cember, the lobes’ central frequencies 
were similar and the large lobe was be- 
tween the radar and the small lobe. No 
ec lips< : of the small lobe was evident, so a 
POS view of rhe asteroid would be expect- 
ed to show both lobes. However, during 10 
to 12 December, when the small lobe was 
closer, the large lobe’s limb extended sev- 
eral microseconds further at frequencies 
not shared by the small lobe than it did at 
frequencies containing echoes from both 
lobes, indicating partial occultation of the 
large lobe by the small lobe. The 1 1 De- 
cember image also shows the echo s great- 
est delay depth. Thus it appears that the 
lobe-lobe line coincided approximately 
with the asteroid’s longest dimension and 
that both were closer to the LOS on 1 1 
December than on any other date. The 
delay dispersions of the two components’ 
echoes place lower bounds of 0.9 and 1.9 
km on their maximum overall dimensions, 
bur further interpretation of the images 
requires knowledge of the kilometers-per- 



Fig. 1 . Toutatis radar images (Table 1). (A) Goldstone low- resolution images and (B) Arecibo images 
plotted with time delay increasing toward the bottom and Doppler frequency increasing toward the left. 
Numbers indicate the date in December 1 992. On the vertical sides, ticks are 2 p.s apart. Two horizontal 
sides have ticks separated by 1 Hz for Goldstone and 0.28 Hz for Arecibo; those intervals correspond to 
a radial velocity difference of 18 mm s Gray scales were designed to show the extent of the 
delay-Doppler distributions. 
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hertz conversion factor as a function of 
time. 

Toutatis’s hour-to-hour rotational mo- 
tion was barely perceptible, but day-to-day 
changes in the orientation of the line con- 
necting the two lobes were obvious. During 
the first few days, the motion of the lobe- 
lobe line and naive application ot Eq. 5 
suggested a rotation period between 4 and 5 
days [compared with 0.1 to 1.0 days for 
almost all other measured ECA rotation 
periods (5)] and a view closer to equatorial 
than pole-on. For example, because the as- 
teroid’s typical overall dimension appears to 
be about 2 km and its typical bandwidth 
was about 2 H: over most of the imaging 
sequence, one might gather from Eq. 5 that 
the POS component ( I W (ipr X e I ) of W. ipp 
was roughly 90° per day (6). However, as 
the experiment progressed, it became in- 
creasingly difficult to identify any constant 
value for W that was consistent with the 


evolution of the delay-Doppler signature, 
and by the end of the experiment w f e real- 
ized that our premise of principal-axis rota- 
tion was incorrect. 

According to force-free rigid-body dy- 
namics, if the angular momentum vector L 
is not aligned with a principal axis, the 
body-fixed inertia ellipsoid will freely pre- 
cess about L while W executes a periodic 
motion in the body, its direction defining a 
closed curve on the inertia ellipsoid (7). Let 
the principal axes x v x 7 , and x ? correspond 
to moments of inertia /j > J 7 > J v Rotation 
is stable only about x ] or x v Convergence of 
L with X| generally will occur whenever 
energy dissipation is taking place in a quasi- 
rigid body (8); light curves show that such 
rotation is predominant among the aster- 
oids (9). In contrast, for Toutatis, the lobes 
always overlap in frequency, even w j hen 
they are at similar ranges, so we never see 
W normal to the longest dimension. Be- 

•>pr G 


cause I W sky I is small compared to I W X 
e I , w'e conclude that W is not normal to the 
asteroid’s longest dimension; in fact, most 
of the images suggest that W is closer to the 
longest dimension than to its normal. This 
circumstance suggests either non-principal 
axis rotation or very unusual physical prop- 
erties, but in the absence of information 
about the asteroid’s shape and internal den- 
sity distribution (and hence the inertia el- 
lipsoid’s shape and body-fixed orientation), 
it does not in itself preclude principal axis 
rotation. 

However, other aspects of the image 
sequence do demonstrate variation of W, 
thereby establishing non-principal axis 
rotation. If W w r ere constant in body-fixed 
and inertial frames, the similarity of the 2 
and 13 December images (for which e 
differs by 110°) would require 8 to be 
similar on each date, but that possibility is 
ruled out by the bandwidth collapse on 1 1 


Table 1 . Toutatis observations, ephemeris corrections, and astrometry. 
Goldstone tmage designations are the day of year followed by a two-digit 
index. Goldstone observations used the 70-m antenna, DSS-14, for trans- 
mission. Negative image designations indicate bistatic observations, done 
by transmitting continuously from DSS- 1 4 and receiving at DSS- 1 3, a 34-m 
antenna 22 km away. Arecibo observations were monostatic. For each 
image in this paper we give Toutatis’s right ascension (RA), declination 
(DEC), distance, and rate of sky motion ( I W sky ! ) at a UTC epoch (f 0 ) near the 
midpoint of data reception; N is the number of looks summed. The frequen- 


cy resolution Af has a radial-velocity equivalent of Af/2 wavelengths per 
second and a length equivalent of Af times the kilometers-per- hertz con- 
version factor from the Hudson-Ostro (74) model. For Arecibo, the data 
processing introduced high correlation between adjacent frequency cells. 
Delay-Doppler ephemeris 27N was calculated on 27 November, 2D on 2 
December, etc. Av eph and A7 eph are model-based (74) COM corrections to 
the specified ephemeris. Astrometry corresponds to echoes received at t 0 
from Toutatis’s COM. For DSS- 13 and DSS-14, the reference point is the 
intersection of the azimuthal and elevation axes. 


Date 

(Dec 

1992) 

^0 

(hh: 

mm) 

Image 

N 

A f 
(Hz) 

km/ 

Hz 

Af x A t 
(m x m) 

RA DEC Dist. 
(deg) (deg) (AU) 

iw sky l 

(°/d) 

Eph. 

^ V eph 

(Hz) 

^ ‘eph 

(m-s) 

Time delay 

(S i (AS) 


Doppler frequency 
(Hz ± Hz) 







Goldstone (8510- 

MHz, 3.5-cm) low- resolution (At = 500 ns) images (Fig. 1) 
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Goldstone (8510-MHz, 3.5-cm) high-resolution (At - 

125 ns) images (Fig. 2) 
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Arecibo (2380-MHz, 

13-cm , At = 200 ns) images (Fig. 1) 
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December, only 13° from the 13 Decem- 
ber position. Also, the 1 1 December band- 
width is about half of that in cither adja- 
cent image, although the asteroid’s direc- 
tion was changing by only 7° from day to 
day. If all other factors were equal and if 
W were constant in body-fixed and iner- 
tial frames, then the twofold change in 
bandwidth within a 7° excursion in e 
would imply that W. ipp was about 7° from 
the LOS on 1 1 December. If we take 2 km 
as an upper bound on the asteroid’s POS 
extent normal to W during these dates, 
then I W ;ipr X e| was at least 17° per day 
on 1 1 December and twice that on the 
adjacent dates. Because l\V skv l was only 
7° per day, the bandwidth collapse on 1 1 
December must have been caused primar- 
ily by the collapse of |W X el, which in 

Table 2. Toutatis’s orbital elements and 3tr un- 
certainties (in parentheses, in units of the last 
decimal places}, estimated with JPL planetary 
ephemeris DE 245 and confirmed with ephem- 
eris PEP740 from the Center for Astrophysics, 
using our radar astrometry (Table 2) and over 
400 optical astrometric observations from 1988 
through April 1993. Angular orbital elements are 
referred to the mean ecliptic and equinox of 
J2000. Weighted rms residuals are 0.83 arc sec, 
0.087 Hz (1 .5 mm s~ 1 in radial velocity) and 0.39 
ixs (59 m in range). TDB, barycentric dynamical 
time. 


Epoch (TDB) 

13.0 Nov 1996 

Perihelion passage time 

11.04640(2) Nov 1996 

(TDB) 

Perihelion distance (AU) 

0.920615832(13) 

Eccentncity 

0.634085830(5) 

Argument of perihelion 

274.45243(135) 

(degrees) 

Longitude of ascending 

128.64567(135) 

node (degrees) 

Inclination (degrees) 

0.46751(1) 


turn could mean either that |W| was 
much smaller than on dates far from 1 1 
December (that is, the magnitude of W 
varies) or that W was very close to the 
lobe-lobe line on 1 1 December. However, 
because the lobe-lobe line itself rotated 
from day to day through most of the rest of 
the sequence, the latter possibility would 
mean that the body-fixed direction of W 
varies. Hence, the image sequence cannot 
he reconciled with a spin vector that is 
constant in inertial and body-fixed frames. 

Burns (10) argued that for plausible as- 
teroid shapes in non-principal axis rota- 
tion, the time scale for precession of the 
inertia ellipsoid about L will be within a 
factor of 10 of the “rotation” time scale (the 
period of W’s body-fixed motion). That 
those time scales are on the order of days for 
Toutatis is suggested by visual inspection of 
Fig. 1 and the above analyses. Because the 
time scale for damping such slow non-prin- 
cipal axis rotation exceeds the age of the 
solar system (9, 1 1 ), Toutatis’s spin state 
might he primordial. Regardless of when it 
was acquired, Toutatis’s rotation offers a 
modern-day, slow-motion example of post- 
collision rotation states likely to have char- 
acterized most asteroids (12). 

Detailed definition of Toutatis’s spin 
state (and shape) requires inversion (13) of 
the image sequence with a comprehensive 
physical model. Hudson and Ostro (14) 
have carried out such an inversion on the 
Goldstone images in Fig. 1 and the 19 
December Arecibo image, and we use some 
of their results here. The model-based 
ephemeris corrections Ai\. ph and Av tph (Ta- 
ble 1) yield astrometry that has been used to 
refine the asteroid’s orbit (Table 2); Touta- 
tis will pass 0.010360 AU (4 lunar distanc- 
es) from Earth on 29 September 2004, the 



Fig. 2. High-resolution Goldstone images (3) of Toutatis from the indicated December dates, 
dilated according to kilometers- per- hertz conversion factors predicted by the Hudson-Ostro 
model (74) for the asteroid’s spin state, with look-to-look translational smearing removed according 
to the model-based estimates of Av eph (f) for each date (Table 1 and Eq. 5). Geometry is as in 
Fig. 1. Perimeter ticks are 1 km apart and gray scales were designed to emphasize brightness 
contrast. 


closest approach predicted for any asteroid 
or comet between now and 2060 (15). The 
corrections also have been used to dc-smear 
high-resolution Goldstone images, which in 
Fig. 2 have been dilated according to kilo- 
meters-per-herti conversion factors from 
the model solution for W ;)pp (r). 

The clearest images in Fig. 2 reveal 
numerous craters with diameters from 
~100 to —600 m; in the delay-Doppler 
projection, concavities are most readily 
discernible near the target’s limb. The sur- 
face density of 100- to 600-m craters on 
Toutatis does not appear to be radically 
lower than that on asteroids 951 Gaspra 
(16) or 243 Ida (17) or on Ida’s satellite 
Dactyl (J8). Inversion of our highest res- 
olution images should define the crater 
distribution in detail. 
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spatial resolution arises from the ability to resolve the echo 
in time delay and Doppler frequency and does not depend 
on the angular size of the target. In practice, such observa- 
tions have been limited by the sensitivity of available tele- 
scopes, and almost all MBA radar measurements are re- 
solved only in Doppler frequency, providing one- 
dimensional images that can be thought of as scans of 
radar brightness taken through a slit parallel to the target’s 
projected spin vector. For the bulk of extant observations, 
the signal-to-noise ratio (SNR) of an optimally filtered sum 
of all spectra obtained in an experiment is less than 10, 
providing marginal spatial resolution. For experiments 
with much higher SNR, the spectra provide unique infor- 
mation about size, shape, and spin vector, but even the 
lowest SNR observations are useful because disc-inte- 
grated radar properties constrain macroscopic roughness 
and near-surface bulk density. 

Ostro, Campbell, and Shapiro (1985, henceforth 
OCS85) reported initial results from Arecibo (13-cm, 
2380-MHz) radar observations of 20 MBAs. Their princi- 
pal results included (i) polarization ratios that indicate, 
for almost all of their targets, a deficit of near-surface 
roughness within an order of magnitude of the wave- 
length, (ii) broad echo spectral shapes that require consid- 
erable roughness at some scale(s) larger than the wave- 
length, and (iii) a fivefold variation in radar albedo, 
implying metal concentrations ranging from nearly zero 
to nearly unity if the surfaces all have regoliths with 
porosities like those on the Moon. Since that paper was 
published, the list of radar-detected MBAs has increased 
by 50% and several objects have been observed during 
more than one apparition. 

Here we present detailed analyses of all radar observa- 
tions of 7 Iris, 9 Metis, 12 Victoria, 216 Kleopatra, and 
654 Zelinda conducted at Arecibo during 1980-1989. 
For Iris, the analysis also incorporates 1991 Goldstone 
radar observations, which were reported by de Pater et 
al (1994). Echoes from each of these asteroids show 
evidence for large-scale topography. For Iris, Metis, and 
Zelinda, Doppler-resolved spectral features appear over 
small intervals of rotation phase, probably due to reflec- 
tions from large flat regions. For Kleopatra and Victoria, 
bimodal spectra in restricted phase intervals indicate 
nonconvex, possibly bifurcated shapes. Our observations 
yield new constraints on pole direction and surface prop- 
erties for all five objects. 

In the next section we describe our observations and 
data reduction and give an overview of our strategies for 
data analysis and physical inference. In Section 3, we pre- 
sent radar results for each asteroid in turn, gearing our 
analysis to the particular combination of available echo 
strength, rotation phase coverage, and nonradar prior in- 
formation about physical properties. In Section 4 we dis- 
cuss the ensemble of results and their implications for 


asteroid science. We conclude by describing immediate 
prospects for observations of these objects with the up- 
graded Arecibo radar. 

2. OBSERVATIONS AND DATA REDUCTION 

Observational and data reduction techniques were very 
similar to those described by Ostro et al (1983 and 1992). 
Briefly, each transmit-receive cycle, or run, yielded echo 
power spectra that we blocked into 4-min sums. Each run’s 
receive interval typically lasted several seconds less than 
the round-trip echo time delay (RTT), between 14 and 
24 min for our observations. Table I lists first and last 
observation dates, total numbers of runs, and average val- 
ues of right ascension, declination, and distance for each 
apparition. The post-1981 observations used a two-channel 
receiving system for simultaneous recording of echoes in 
the same sense of circular polarization as transmitted (the 
SC sense) and in the opposite (OC) sense. (The 1980 Iris 
observations used a single-channel system, which was 
switched between OC and SC in alternate runs.) The heiic- 
ity of circular polarization is reversed upon reflection from 
a surface that is smooth on all scales within about an order 
of magnitude of the wavelength, but SC echo power can 
arise from single backscattering from a rough surface, from 
multiple scattering, or from subsurface refraction. The cir- 
cular polarization ratio /x c = SC/OC is thus a measure of 
the near-surface structural complexity, or ‘Toughness,” at 
scales near the wavelength. 

Echo power is given by P R = P T G T G R A 2 cr/(47r) 3 /? 4 , 
where P r is the transmitted power, G T and G R are the 
antenna gains during transmission and reception, A is the 
observing wavelength, R is the radar-target distance, and 
cr is the radar cross section (either OC or SC), defined as 
4rr times the backscattered power per steradian per unit 
incident flux at the target. Echo power was measured as 
a function of frequency relative to the Doppler frequency 
of hypothetical echoes from the center of mass (COM) as 
predicted by site ephemerides. The uncertainties in the 
Doppler-prediction ephemerides were small compared 
with the data’s frequency resolution. 

We normalized each of our 4-min spectra to the standard 
deviation of the receiver noise to facilitate the formation 
of weighted sums of spectra taken with different values 
of antenna gain, transmitter power, system temperature, 
target distance, and integration time. Uncertainties in esti- 
mates of radar cross section are due primarily to systematic 
errors in calibration of antenna gain (as a function of eleva- 
tion and azimuth), transmitter power, and system tempera- 
ture, as well as pointing accuracy, which varies over a 
variety of time scales. Our experience during the 1980s 
with observations of a variety of radar targets leads us to 
believe that for the most part, the absolute uncertainty in 
radar cross section estimates is between 20 and 50%, while 
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TABLE I 
Observations 




to 

No. 

RA 

Dec 

Dist 

/tx 

A/ 

Target 

Dates spanned (UT) 

(hr. UT) 

runs 

(hr) 

(deg) 

(AU) 

(MHz) 

(Hz) 


1980 Sep 26-Oct 01 

2.62 

7 

23.6 (0.07) 

10 (0.6) 

0.94 (0.001) 

2380 

9.8 

7 Iris 

1984 Dec 19-Dec 24 

2.84 

25 

5.3 (0.08) 

22 (0.5) 

0.95 (0.018) 

2380 

9.8 


1991 Sep 15-Sep 28 

4.43 

32 

22.7 (0.17) 

3(1) 

1.02 (0.012) 

8495 

39.1 

9 Metis 

1984 Mar 17-Mar 22 

3.39 

15 

12.2 (0.08) 

9 (0.4) 

1.44 (0.006) 

2380 

9.8 

1986 Nov 21 -Nov 23 

2.25 

9 

4.4 (0.04) 

21 (0.04) 

1.12 (0.003) 

2380 

9.8 

12 Victoria 

1982 Sep 29-Oct 05 

1.16 

15 

0.2 (0.08) 

15 (1.1) 

1.03 (0.014) 

2380 

9.8 

1989 Aug 19-Aug 21 

8.62 

9 

22.3 (0.03) 

9 (0.1) 

0.88 (0.001) 

2380 

9.8 

216 Kleopatra 

1985 Nov 18-Dec 01 

0.21 

13 

4.0 (0.16) 

9 (2.0) 

1.18 (0.024) 

2380 

18.8 

654 Zelinda 

1988 Jan 17-Jan 18 

6.70 

9 

5.7 (0.01) 

21 (0.3) 

0.89 (0.004) 

2380 

3.9 


Note. Right ascension, declination, and distance are given for epochs near the weighted midpoint of observation; the range of values spanned is 
in parentheses. An epoch of zero rotation phase, / 0 , is given for the first date of each experiment. Zero phase does not correspond to any particular 
orientation of the asteroid, nor is there any phase correspondence from year to year. The number of transmit-receive cycles, or runs, is given in 
the fourth column. The transmitter frequency (/ TX ) and the raw frequency resolution (A/) of the echo Doppler spectra are given in the last two columns. 


relative uncertainties are half as large. Most systematic 
effects cancel for fi c > and statistical uncertainty from the 
propagation of receiver noise (Appendix I of Ostro et al 
1983) dominates our quoted errors. The time-bandwidth 
products for all spectra presented in this paper are >100, 
and the noise obeys Gaussian statistics. Our method for 
estimation and removal of the mean background of the 
receiver noise was described most recently by Ostro et al. 
(1992). Typically, several runs were obtained each day 
over a span of several days to provide spectra at different 
rotation phases and to increase the total SNR. For each 
radar apparition, the target’s plane-of-sky motion was very 
small (Table I). The useful rotation phase resolution is 
restricted by the available sampling in rotation phase and 
the echo SNRs and is different for each target. 

3. DATA ANALYSIS 

An echo’s strength and polarization and its distribution 
in Doppler frequency and time delay contain information 
about the target’s characteristics on size scales from its 
overall dimensions down to the microscopic. It is useful 
to imagine a hierarchy of scales 

wavelength facet tile < asteroid, 

where “tile” implies the effective spatial resolution of the 
data and/or an assumed model. Most radar-detected aster- 
oids, including those discussed in this paper, have low circu- 
lar polarization ratios (/ x c ~ 0.1), which indicates that most 
of the OC echo power is due to single back-reflections 
from surface elements, or “facets,” that are smooth on all 


scales within about an order of magnitude of the wave- 
length (OCS85). 

For targets with low /i C , the OC radar albedo (doc = 
a- 0C Mp, where A v is the target’s projected area) can be 
related to the Fresnel power reflection coefficient at normal 
incidence ( R ) via: 

Ooc - = gR, (1) 

where the gain factor g depends on the target’s orientation, 
its gross shape, the distribution of surface slopes with re- 
spect to that shape, and the degree of wavelength-scale, 
near-surface roughness. For a sphere with /x c = 0, the gain 
factor is unity, and one can think of g's departure from 
unity as quantifying how the distribution of facet tilts with 
incidence angle differs from that of a sphere. For most 
large main-belt asteroids, g is expected to be within a few 
tens of percent of unity (OCS85), so that a oc provides a 
reasonable first approximation to R. For dry, particulate 
mixtures of rock and metal with particle sizes <A/100, R 
depends strongly on bulk density, which in turn is a func- 
tion of the porosity, the metal weight fraction, and the 
specific gravities of the rock and metal phases. 

An echo’s instantaneous edge-to-edge bandwidth B at 
rotation phase 4> can be written 

B(<M) = ^§^cos(<S), (2) 

where D is the breadth normal to the radar line of sight 
of the asteroid’s pole-on silhouette, S is the target-centered 
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declination of the radar, and P is the apparent rotation 
period. Each target's rotation period is known quite well 
(Lagerkvist et al 1989 and references therein) and the 
contribution of plane-of-sky motion to the apparent rota- 
tion was negligible for all of our observations; therefore, 
we treat our bandwidth estimates as joint constraints on 
D(<f>) and S. Our ability to discern the spectral edges of an 
echo and hence the perceptibility of rotational bandwidth 
variations depends on the shape and radar scattering prop- 
erties of the target and on the SNR, which for the targets 
discussed here is not very high. In any event, the observed 
bandwidth can never be greater than the maximum band- 
width, B max (S) = [47rD max / \P] cos (£), corresponding to 
the maximum breadth of the target's pole-on silhouette. 

One goal of asteroid astronomy is the development 
of realistic, detailed shape models. For several of our 
targets, ellipsoid models based on photometric 
lightcurves, occultations, and infrared radiometry provide 
approximations that refine radar albedo estimates and 
also provide a priori estimates of 5 max (<5). When there 
is sufficient confidence in knowledge of an echo’s band- 
width, the spectral shape constrains surface structural 
characteristics at both facet and tile scales. For example, 
one measure of surface roughness on scales 5:10A is the 
fractional half-power bandwidth B HP IB m ^ x . Radar echoes 
from the Moon, with B HP iB max — 0,1, are dominated by 
specular reflections from surface elements near the sub- 
radar point. OCS85 found that B HP !B max — 0.5 for the 
very-low -p c targets Ceres and Pallas, indicating surfaces 
that are much rougher than the lunar surface at some 
scale(s) greater than a few meters. 

With adequate SNR, echo bandwidth and spectral 
structure as functions of rotation phase can constrain an 
asteroid’s shape by revealing topographic features. Higher 
SNR cases may warrant using the echoes to estimate 
the shape. The number of shape parameters constrained 
by the data is dictated by the SNR, the Doppler resolu- 
tion, and the orientation coverage of the data. Some 
approaches used to model Earth-crossing asteroid echoes 
(e.g., Ostro et al 1990, Hudson and Ostro 1994) have 
limited applicability to the observations reported here. 
All such approaches parameterize the surface’s average 
angular scattering law [a- o (0) = daldA , where dA is an 
element of surface area and 6 is the angle of incidence], 
which quantifies effects of structure at facet scales. Simp- 
son and Tyler (1982) described efforts to infer facet 
slope statistics from measurements of a o (0) within the 
framework of various definitions of slope probability 
density functions. For example, in modeling echoes from 
Pallas, OCS85 assumed a spherical shape and a scattering 
law proportional to exp(-So 2 tan 2 (0))/cos 4 0, where s 0 
is the adirectional rms slope, Parker’s (1973) one-dimen- 
sional probability density function is implicit, and 
Gaussian height distribution and surface autocorrelation 


functions are assumed. All other modeling of asteroid 
echoes to date has used an empirical cos" (6) law, for 
which n = 1 corresponds to geometric scattering (radar 
brightness proportional to projected area), n ~ 2 to 
Lambert limb darkening, and higher values to more 
specular backscattering. Estimates of the radar scattering 
law exponent for two Earth-crossing asteroids yield 
n ^ 2 for 1685 Toro (Ostro et al 1983) and n - 
2.8 ± 0.3 for 4769 Castalia (Hudson and Ostro 1994). 
If a Parker probability density function is assumed, 
then a cos" (0) scattering law implies that s 0 = \ // 2/n 
or tan _1 (s 0 ) = 45 and 40° for those two objects (see the 
Appendix). 

4. RESULTS 

In this section, we begin the discussion of each target 
by assessing the available nonradar information about size, 
shape, and spin vector from radiometric, lightcurve, and 
occultation observations. These constraints are summa- 
rized in Table II for ease of reference. For Iris, Metis, and 
Kleopatra, analyses of multi-apparition lightcurves have 
produced estimates of pole direction as well as simple 
shape constraints in the form of model ellipsoid axis ratios 
(a/b and b/c , where a > b ^ c and rotation is about the 
shortest axis). Size constraints derived from a combination 
of optical and infrared data are available for Iris, Victoria, 
Zelinda, and Kleopatra (Tedesco and Veeder 1992). Occul- 
tation observations contribute prior information about the 
sizes and shapes of Metis and Kleopatra. 

Galileo’s encounters with Gaspra and Ida have provided 
two tests for evaluating groundbased predictions for size, 
shape, and spin vector. Although both asteroids are irregu- 
larly shaped, ellipsoid models derived from groundbased 
VIS/IR observations furnished accurate predictions of size, 
elongation, and (ambiguities aside) pole direction (Mag- 
nusson et al 1992, Binzel et al 1993, Belton et al 1994). 
We have studied the published nonradar constraints on 
our targets’ geometric properties and in several cases have 
settled on an a priori ellipsoid model as the starting point 
of our analysis. 

7 Iris 

Magnusson’s (1986) analysis of Iris lightcurves yielded 
a model ellipsoid with a/b - 1.18, b/c = 1.41, and a pole 
direction with ecliptic longitude and latitude (A, fi) of either 
(15 ± 5°, +25 ± 15°) or (195 ± 5°, +15 ± 15°). The 
lightcurve coverage for Iris is such that uncertainties in 
the axis ratios are believed to be <0.1 (Magnusson 1990). 
Another lightcurve analysis (Zappala and Di Martino 
1986) yielded a model ellipsoid with a/b — 1.19, b/c = 
1.21, and possible pole directions of (18 ± 5°, +33 ± 7°) 
and (193 ± 4°, +16 ± 8°). These analyses are in agreement, 
except for their estimates of b/c , possibly because of differ- 
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ent assumptions about the optical scattering properties 
of Iris’ surface. Magnusson assumed geometric scattering, 
while Zappala and Di Martino applied an empirical correc- 
tion factor based on the laboratory work of Barucci et al. 
(1984). Magnusson (1990) argues convincingly that geo- 
metric scattering is a better approximation to the true 
scattering law than is the model of Barucci et al., so we 
will adopt Magnusson’s model. 

Goldstone-VLA radar aperture synthesis observations 
of Iris (de Pater et al. 1994) removed the ambiguity in the 
pole direction in favor of the solution (A = 15°. B = +25°). 
This pole direction corresponds to a subradar latitude 
S = 19 ± 8° during the IRAS observations, which provided 
a radiometric diameter of 200 ± 10 km (Tedesco and 
Veeder 1992). The average projected area of Magnusson’s 
ellipsoid model during those observations would be the 
same as that of a 200-km sphere if the ellipsoid’s dimen- 
sions were 260 X 220 X 155 km. We adopt this as a working 
model for Iris, with a 15% uncertainty in each dimension, 
which is intended to include the quoted uncertainties in 
the radiometric diameter and the model ellipsoid s axis 
ratios, possible systematic bias in scaling the ellipsoid using 
the radiometric diameter (Brown 1985), and uncertainty 
in the asteroid’s orientation at the time of the IRAS obser- 

vations. . 

Figure 1 shows weighted sums of echo spectra obtained 
in 1980 and 1984 smoothed to a resolution of 30 Hz, Fig. 

2 reproduces a weighted-sum spectrum from recently re- 
ported 1991 Goldstone 3.5-cm observations (de Pater et 
al. 1994) smoothed to 240 Hz, and Figs. 3 and 4 show sums 
of spectra within selected rotation phase groups for 1980 
and 1984. As depicted in polar plots inset into the figures, 
each spectrum represents a weighted average of spectra 
from short integration times (4-min blocks or individual 
runs — see figure legends). In the polar plots, the standard 
deviation of the noise for each spectrum that contributes 
to the sum is represented by a radial error bar at the 
corresponding relative rotation phase. The rotation phase 
origins are arbitrary and do not correspond to any particu- 
lar orientation of the model ellipsoid. The SNRs of the 
spectra in Figs. 3 and 4 preclude detection of rotational 
bandwidth variations as small as those predicted by the 
model ellipsoid. 

The horizontal extents of the shaded boxes in Figs. 1 and 
2 encompass interval estimates for spectra edge locations 
based on B max (S) for the nominal ellipsoid model and pole 
direction (Table II) and the uncertainties therein. The ap- 
parent extents of the OC spectra (i.e., the bandwidths 
between first zero crossings) are consistent with the prior 
constraints. These spectra, therefore, are best described 
as very broad (B HP IB > 0.5), and we can conclude with 
reasonable confidence that Iris is rough at some scale(s) 
no smaller than the wavelength. However, Iris’ low circular 
polarization ratio (Table III) indicates that most of the 


echo power arises from single back-reflections from facets 
whose sizes and radii of curvature are large compared to 
A. Thus the detection of echo power near the spectral 
edges, far from the center of the asteroid’s disc, seems to 
require considerable roughness at some scale(s) greater 
than a few meters. 

Our 1984 Iris spectra show evidence for roughness at 
topographic scales, that is, scales that are not negligible 
compared with the asteroid’s overall dimensions. A narrow 
spike centered near —300 Hz can be seen in the group 
1” spectrum of Fig. 4. The available SNR does not allow 
precise estimation of the attributes of this spike, but we 
have determined the optimal filtering in rotation phase 
and Doppler frequency that maximizes the spike’s SNR. 
An average of six consecutive spectra spanning a rotation 
phase interval of A<f> = 13°, smoothed to 90 Hz resolution, 
yields SNR « 8 for the spike. Because of limited SNR and 
a data gap immediately following this interval, we consider 
this estimate to be only a crude one for the feature’s extent 
in rotation phase. Figure 5 shows this spectrum and simi- 
larly processed spectra before and after the appearance of 
the radar spike. The 90-Hz frequency filter corresponds to 
a distance, measured normal to the plane containing Iris 
spin vector and the radar line of sight, of ~25 km, which 
sets a lower bound on one dimension of the region respon- 
sible for the spike. If this feature were due to a region of 
unusually high intrinsic radar reflectivity, we would expect 
to see it migrate from left to right as rotation brought it 
from positive to negative Doppler frequencies. The fea- 
ture’s presence in a restricted phase interval indicates that 
the surge in brightness is due to geometry, i.e., a temporary 
surge in surface area oriented normal to the radar. That 
is, there must be a large region on Iris with surface facets 
more or less parallel to each other. Various detailed con- 
figurations of facets are obviously possible, but the most 
plausible physical picture is of a continuous, fairly flat 
region. 

Table III shows estimates of Iris’ OC radar cross section, 
OC radar albedo, and circular polarization ratio from the 
1980, 1984, and 1991 experiments. The radar cross section 
(opc) is calculated from a weighted average of available 
spectra, and the radar albedo (ooc) is estimated by dividing 
ctoc by the projected area of the nominal a priori ellipsoid 
averaged over all rotation phases. The quoted uncertainties 
(in parentheses) include all known sources of error, includ- 
ing the fact that the phase for the ellipsoid’s minimum- 
bandwidth orientation is not known. The weighted average 
of all radar albedo estimates for Iris is cr oc = 0.11 ± 0.03. 
This and values for our other four targets will be discussed 
in Section 4. The 13-cm estimate of Me from 1980 is less 
than either the 3.5-cm estimate at a comparable subradar 
latitude or the 1984 estimate obtained with a more equato- 
rial view. Discrepancies in estimates of mc for other multi- 
year asteroid radar data sets are rare, even at the low level 
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FIG. 1. Arecibo radar spectra of Iris obtained in 1980 (top) and 1984 (bottom). Echo power, in standard deviations of the noise, is plotted 
versus Doppler frequency relative to that of hypothetical echoes from the asteroid’s center of mass. The solid and dotted lines plot the OC and 
SC echoes, respectively, smoothed to a frequency resolution of 30 Hz. Each spectrum is a weighted sum of independent spectra obtained at different 
asteroid rotation phases, which are depicted in the inset with a radial "error bar" proportional to the standard deviation of each spectrum included 
in the average. The arrow indicates zero phase, as defined in Table I, and phase increases in the counter-clockwise sense. There is no correspondence 
between 1980 and 1984 phases. The vertical extent of the shaded boxes shows standard deviation of the noise, while the horizontal extent shows 
the expected range for the edge-to-edge bandwidth (B^) based on radiometric and Iightcurve constraints on the asteroid’s size, shape .and spin 
vector (see Table II). 



Doppler Frequency (Hz) 


FIG. 2. Goldstone 3.5-cm radar spectrum of Iris obtained in 1991, 
smoothed to a frequency resolution of 240 Hz. (See legend to Fig. 1.) 


of significance here, so we take these numbers to suggest 
the possibility of either regional or scale-dependent varia- 
tions in small-scale structure. 

9 Metis 

Since the mid 1980s, three separate Iightcurve analyses 
(all assuming geometric optical scattering) have provided 
triaxial ellipsoid models and pole directions for Metis. Zap- 
pala and Knezevic (1984) found axis ratios of a/b = 1.32 
and b/c = 1.34 and possible pole directions with ecliptic 
coordinates of Z1 = (2 ± 6°, +26 ± 6°) and Z2 = (186 ± 
8°, +43 ± 6°). Magnusson (1990) found axis ratios of 
a/b = 1.27 and b/c = 1.26 and a pole direction of either 
Ml = (0 ± 5°, +20 ± 10°) or M2 = (180 ± 5°, +30 ± 
10°). As for Iris, uncertainties in Magnusson’s axis ratios 
are believed to be <0.1. Finally, Drummond et al (1991) 
obtained axis ratios of a/b = 1.27 ± 0.02 and b/c = 
1.24 ± 0.02 and a pole direction within 8° of either D1 = 
(0°, +7°) or D2 = (181°, +23°). The model axis ratios are 
within Magnusson’s quoted uncertainties, so we will adopt 
Magnusson’s axis ratios with conservative uncertainties of 
0.1. The pole longitudes agree to well within the quoted 
uncertainties, but there is a significant dispersion in the 
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TABLE n 
Prior Information 


Target 

Class® 

IRAS 

diameter* 

Model ellipsoid^ 

Synodic 

period 1 * 

Pole direction* 

Year' 

<5* 

(Ap)* 

B m „ (sy 

2a 

2b 

2c 

long. 

lat. 

7 Iris 

S 

200 ± 10 

260 

220 

155 

7.139 

15 * 5 

+25 ± 15 

1980 

-69 ± 10 

43100 ± 9300 

360 * 180 










1984 

-21 ± 8 

31900 ± 7700 

940 * 150 










1991 

-56 ± 10 

40900 * 9100 

2010 * 620 

9 Metis 

s 

— 

215 

170 

135 

5.078 

0*5 

+20 ± 10 

1984 

+61 ± 10 

27000 * 5800 

570 * 200 








180 + 5 

+ 30 ± 10 


-69 * 11 

27600 * 6000 

420 ± 220 










1986 

-20 * 7 

21500 ± 5400 

1100 * 170 











+20* 7 

21500 ± 5400 

1100 - 170 

12 Victoria 

s 

113 ± 3 

132 

105 

105 

8.662 

9 ± 9 

+55 ± 10 

1982 

+47 ± 10 

10300 * 2300 

290 * 100 








176 ± 4 

+40* 5 


-37 * 6 

10200 * 2300 

340 + 80 










1989 

+47 * 11 

10300 * 2300 

290 ± 100 











-30 * 6 

10100 i 2300 

370 ± 80 

216 Kleopatra 

M 

135 i 2 

250 

90 

70 

5.385 

71 ± 3 

+ 19 ± 3 

1985 

-58 ± 4 

15900 * 3500 

680 * 130 

654 Zelinda 

c 

127 i 4 

— 

— 

— 

31.9 

- 

— 

1988 

- 

12700 ± 2500 

<140 


® Taxonomic classification (Tholen 1989) based on visual and infrared data. 

*Tedesco and Veeder (1992). These “radiometric diameters” (in km) assume a spherical asteroid that obeys the Standard Thermal Model (Lebofsky and Spencer 
1989). Departures from sphericity may bias these estimates to an extent that depends on the asteroid’s shape and its orientation at the time of the IRAS observation 
(Brown 1985). 

c Adopted axis dimensions (in km) based on a combination of radiometric, lightcurve, and occultation data (see text). We adopt an uncertainty of 15% in each dimension. 

d Except for Victoria (Erikson 1990), synodic periods (in hours) are from Lagerkvist et al. (1989) and references therein. See also Lagerkvist et al. (1987, 1988). 

r Ecliptic coordinates (in degrees) of the spin vector. Except for Victoria (Dotto et at 1995), all pole directions are those of Magnusson (1986, 1990). Alternative pole 
directions, which may differ slightly from those listed, are noted in the text. 

f Year of radar observation. All observations were performed with the 12.6-cm Arecibo radar, except for the 1991 Iris observations (de Pater et al. 1994),which were 
performed with the 3.5-cm Goldstone radar. 

K Subradar latitude over the duration of observation corresponding to the listed pole direction and its uncertainties. 

h Unweighted average projected area (in km 2 ) of the model ellipsoid over all rotation phases for the duration of observation. The uncertainty in the model’s projected 
area incorporates the uncertainties in the axis dimensions and pole direction. 

' Maximum edge-to-edge bandwidth (in Hz) over the duration of observation for the a priori ellipsoid with the listed pole direction and synodic spin period and al! 
associated uncertainties. There is no a priori ellipsoid for Zelinda, so we assume the maximum breadth sampled is between 1.0 and 1.2 times the IRAS diameter. 


pole latitudes, which may be due in part to solution re- 
finements as additional lightcurve data became available. 

Metis was not observed by IRAS, but occultation chords 
from 1984 are fit by an ellipse with overall dimensions of 
210 x 170 km (Kristensen 1984; see also Millis and Dun- 
ham 1989). The sky projections of the poles listed above 
are aligned within 40° of the occultation ellipse’s long axis, 
suggesting some foreshortening of the asteroid’s longest 
dimension. Given the above pole directions and their un- 
certainties, the view during the occultation was sufficiently 
close to pole-on (49° < |<5| < 83°) that the projection of 
Magnusson’s ellipsoid would be at least 93% as long as the 
ellipsoid itself. In light of this information, we adopt as an 
a priori shape model an ellipsoid with dimensions within 
15% of 215 x 170 x 135 km. 

Figure 6 shows weighted sums of echo spectra obtained 
in 1984 and 1986. Predicted ranges for fl max based on the 
adopted a priori shape model and the pole directions listed 
above are shown by shaded boxes for comparison with the 
apparent spectral edges (the innermost zero crossings) for 
each year. The various pole solutions give identical predic- 
tions for S and hence B max for the 1986 apparition, but 
different predictions for 1984. Whereas the 1984 spec- 
trum’s edges seem to favor the Dl, D2, and M2 pole solu- 


tions, an unusual shape and/or a highly specular radar 
scattering law could conceivably produce a spectral shape 
with weak wings that would be obscured by the noise. In 
any event, the currently available radar data cannot resolve 
the pole ambiguity for Metis. 

Figures 7 and 8 show sums of spectra within selected 
rotation-phase groups for each year. As with Iris, a radar 
spike in Metis spectra from the 1984 aparition is seen within 
a limited rotation phase interval. This spike dominates the 
“group 5” spectrum in Fig. 7 and is sufficiently strong that 
it can be detected in two adjacent rotation phase subgroups 
(Figs. 9a and 9b) spanning —24° of rotation phase. How- 
ever, the SNR is insufficient to detect the spike at rotation 
phases immediately adjacent to that range (e.g.. Fig. 9c), 
so we cannot discern rotation phase boundaries for the 
feature. The spike’s height is maximized by filtering to a 
resolution of 50 Hz, which corresponds to a linear dimen- 
sion of at least 9.2 km/cos S. The M2 pole direction predicts 
S ~ -69° ± 11°, which corresponds to 26 (+27, -9) km 
or 12 (+27, -9) % of the a priori ellipsoid’s longest dimen- 
sion (vs —25 km and —10% for the Iris feature). The Dl 
and D2 pole directions predict \S\ = 75 ± 11°, which corre- 
sponds to 36 (+95, —15) km or 17 (+44, -7) % of the a 
priori ellipsoid’s longest dimension. Our estimates of /x c 
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TABLE III 

Radar Properties by Experiment 


Target 

Year 

OC SNR' 2 

B eq (Hz)* 

BiipiBz c c 

<7bc (km 2 )'' 

Me' 

&oc? 

Iris (A) 
Iris (G) 

1980 

1984 

1991 

22 

28 

19 

277 

643 

1476 

0.70 

0.54 

0.55 

5900 (1500) 
2900 (700) 
4100 (1000) 

0.08 (0.03) 
0.19 (0.03) 
0.19 (0.05) 

0.14 (0.05) 
0.09 (0.03) 
0.10 (0.03) 

Metis 

1984 

1986 

12 

13 

180 

640 

0.86 

0.61 

3500 (920) 
2900 (730) 

0.18 (0.08) 
0.13 (0.05) 

0.13 (0.04) 
0.13 (0.05) 

Victoria 

1982 

1989 

33 

54 

180 

160 

0.56 

0.60 

2100 (520) 
2500 (620) 

0.14 (0.03) 
0.13 (0.05) 

0.20 (0.07) 
0.24 (0.08) 

Kleopatra 

1985 

19 

570 

0.69 

7100 (1800) 

0.00 (0.05) 

0.44 (0.15) 

Zelinda 

1988 

37 

80 

0.56 

2300 (580) 

0.13 (0.03) 

0.18 (0.06) 


“The OC SNR is the signal-to-noise ratio for an optimally filtered, weighted sum of all spectra from an experiment, 

o Bv definition (Tiuri 1964 ), B EO = 4jF[<2«*/2Sfl, where S, are the spectral elements and A/ is the “raw” frequency resolution. 

c sir and B zc are the half-power and zero-crossing bandwidths, respectively, of the weighted sum of all spectra from an experiment smoothed 

As,i e „,d «. ,h. too. mm K,U,„ of systematic ciibt.ti.o e.ot, estimated as 25* of the 

cross-section values, and the standard deviation of the receiver noise in the equivalent bandwidth (Be o)- 

« „ c is the circular polarization ratio (of SC to OC echo power). The standard deviations quoted for mc propagate from the receiver noise . 
/The radar albedo froo is obtained by dividing ^ by the average projected area of the <r priori model ellipsoid at the epoch of our observations, 
if available or the projected area of a sphere with the IRAS diameter (see Table II). Uncertainties propagate from those given for croc and <A P >. 





+600 0 -600 


FIG S OC spectra of Iris obtained in 1984 within narrow rotation phase ranges (cf. Fig. 4). A radar spike appears in spectrum "b centered 
at a Doppler frequency of -305 Hz. but not in adjacent rotation phase ranges (a and c). The spectra are smoothed to a frequency resolution of 
Hz to maximize the spike’s height. The shaded boxes are as in the lower panel of Fig. 1. 
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FIG. 6. Arecibo radar spectra of Metis obtained in 1984 (top) and 1986 (bottom). The solid and dotted lines plot the OC and SC echoes, 
respectively, smoothed to frequency resolutions of 50 Hz (top) and 150 Hz (bottom). The notation is the same as in Fig. 1. Various lightcurve- 
based pole situations (see text) result in different predictions for the edge-to-edge bandwidth in 1984 (stacked shaded boxes), but identical predictions 
in 1986. The innermost zero crossings of the 1984 spectrum seem to favor the Dl, D2, and M2 pole solutions (see text). 


for Metis (Table III) are similar to those for Iris, so our 
inferences about surface structure on Iris seem broadly 
applicable to Metis. However, the Metis glint was seen 
during a more pole-on apparition (|£| Metls - 50°) and the 
Iris glint was seen during a more equatorial apparition 
(<5jri s = 19 ± 8°). Therefore, the putative flat regions respon- 
sible for the glints must have very different orientations 
with respect to the two asteroids’ equatorial planes. 

654 Zelinda 

IRAS observations of Zelinda yielded a radiometric di- 
ameter of 127 ± 4 km, and lightcurves have established a 
rotation period of 31.9 hr but no a priori shape models 
or pole-direction constraints have been published for this 
object. Nonetheless, Zelinda’s lightcurves exhibit peak-to- 
valley brightness variations of —0.3 magnitudes (Schober 
1975). Noting that Iris’ lightcurves have similar amplitudes 
and that projected area variations of Iris’ a priori ellipsoid 
model are as large as 18%, we adopt a 20% uncertainty 
for Zelinda’s projected area. 

We observed Zelinda on Jan. 17 and 18, 1988, completing 
four runs on the first date and five on the second. Each 
run yielded four OC/SC pairs, one from each of four 4- 
min accumulations. Figure 10 plots OC weighted-sum spec- 
tra for (a) the nine runs, (b) the two dates, and (c) the 
entire experiment. The two dates spanned two sections of 


rotation phase about 82° apart; 166-191° on Jan. 17 and 
93-100° on Jan. 18. Our phase origin is arbitrary — we do 
not know the correspondence between our phases and 
those of lightcurve extrema. 

For a target with Zelinda’s rotation period, Eq. 2 can 
be rewritten to relate instantaneous bandwidth ( B , Hz) to 



FIG. 7. Weighted sums of OC echo spectra of Metis from 1984 within 
the rotation phase intervals depicted in the polar plot, each smoothed 
to a frequency resolution of 60 Hz. The frequency scale below spectrum 
6 applies to all spectra. 
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FIG. 8. Weighted sums of OC echo spectra of Metis from 1986 within 
the rotation phase intervals depicted in the polar plot, each smoothed 
to a frequency resolution of 100 Hz. The frequency scale below spectrum 
2 applies to all spectra. 


breadth (D, km): B = 0.869 D cos 8 , so estimation of B 
places a joint constraint on D and 8. The innermost zero 
crossings of the spectrum in Fig. 10c give B = 98 Hz, 
corresponding to D cos 8 = 113 km. A more sophisticated 


estimator uses least squares to fit a model spectrum to the 
data, i.e., to extrapolate from echoes above the noise down 
to the spectral edges. A very simple model, corresponding 
to a spherical target with a cos" 0 radar scattering law, has 
the form: 



Here, the spectral shape is determined by a single param- 
eter, ft, so that for nonspherical targets, the effects of target 
shape, orientation, and scattering properties are all ab- 
sorbed by n. A correlation between B and n is “built in” 
to this model: echo spectra of larger, more specular targets 
may be difficult to distinguish from those of smaller, less 
specular ones. 

Fits of an S(f) model yield n = 2.4 ± 0.7 and B = 
89 ± 4 Hz. This value of B, which corresponds to D cos 
8 = 102 ± 5 km, is indicated in Fig. 10c, along with B(n) 
for 1 < n < 6, a range that spans values of n reported so 
far for fits of an S(f) model to asteroid radar data, and 
also for modeling in which the target shape is so thoroughly 
parameterized that one can realistically interpret n as a 
measure of scattering specularity. For example, in their 
analysis of echo spectra of the rather elongated asteroid 



+500 0 -500 



+500 0 -500 +500 0 -500 


FIG. 9. OC spectra of Metis obtained in 1984 within narrow rotation phase intervals (cf. Fig. 7). A radar spike appears in spectra “a" and “b" 
centered at a Doppler frequency of —85 Hz. The spectra are smoothed to a frequency resolution of 50 Hz in order to maximize the spike’s height. 
No statistically significant feature is present in spectrum “c,” which has a higher noise level than the other two spectra. The shaded boxes show 
edge-to-edge bandwidth predictions based on the D1 and D2 pole solutions (see Fig. 6). 
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FIG. 10. Radar Doppler spectra of Zelinda obtained at Arecibo in 1988. (a) Weighted average OC spectra from each run (four 4-min 
accumulations) smoothed to a frequency resolution of 12 Hz. The weighted average rotation phase is indicated above each plot, (b) Weighted 
average OC spectra for each day at the raw 4-Hz frequency resolution, (c) Weighted average OC and SC spectra for the entire experiment at the 
raw 4-Hz frequency resolution. The horizontal extent of the shaded boxes is derived from fitting the spectrum with an S(f) model (Eq. 3) for 1 < 
n — 6* The best-fit value of n is 2.4 ± 0.7 (arrows labeled “L.S.”). The arrows labeled “IRAS” show' the maximum edge-to-edge bandwidth 
corresponding to the IRAS diameter and an equatorial view. 


1685 Toro, Ostro et al. (1983) obtained n = 3.6 ± 1.4 from 
fitting an S(f ) model spectrum to a weighted sum of 47 
spectra with thorough rotation phase coverage, and n = 
2.04 ± 0.45 from fitting the entire data set with a biaxial 
ellipsoid model with a cos" 0 scattering law. Also indicated 
in Fig. 10c is the bandwidth, 110 Hz, for an equatorial view 
of a sphere with Zelinda’s IRAS diameter. Given Zelinda’s 
lightcurve amplitude (0.3 mag), it seems plausible that the 
asteroid’s maximum breadth might be —15% larger than 


the IRAS diameter. We don’t know if our data sample 
the maximum-breadth orientation, but surely we sampled 
orientations <45° from it. In consideration of all this infor- 
mation, we conclude that our view of Zelinda was probably 
closer to equatorial than to pole-on. 

The Jan. 18 spectra exhibit a brightness feature that 
persists from 92.4 through 98.8° (Fig. 11), but because of 
a gap in phase coverage from 85 to 92°, this feature could 
be present at a similar strength for up to 14° of rotation 
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FIG. 11. OC spectra of Zelinda from individual 4-min accumulations obtained at Arecibo on 18 January 1988 (cf. Fig. 10). The rotation phase 
is indicated above each spectrum. All spectra are smoothed to a frequency resolution of 12 Hz, and the frequency scale below the last specrtrum 
applies to all spectra. A radar spike appears at 92.4° and disappears abruptly at 99°. 
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phase. The bandwidth of the feature is —12 Hz (three 
spectral resolution elements), corresponding to a linear 
dimension of at least 14 km/cos S. As with the Iris and 
Metis glints, the source cannot be a reflectivity feature but 
rather must arise from an anomalously large fraction of 
the asteroid’s projected area in a small Doppler domain 
being normally oriented toward the radar. Even A <f> — 7° 
is much larger than the <^1° range that would result from 
a perfectly flat, kilometer-sized plane rotating through the 
radar-facing orientation. The observed phase-frequency 
signature of the radar feature probably could result from 
a region whose surface elements have normals within -10° 
of each other. We interpret the rapid fading of the glint 
at 99° to mean that the responsible region is extremely flat 
at scales very much larger than 10 m. Zelinda’s low circular 
polarization ratio indicates that the surface is very smooth 
at centimeter-to-meter scales. 

216 Kleopatra 

Kleopatra’s maximum optical lightcurve amplitude of 
1.18 magnitudes is among the largest measured for any 
main-belt object (Lagerkvist etal 1989) and is remarkable 
considering the asteroid’s large radiometric diameter, 
135 ± 2 km. Extensive lightcurve observations have been 
interpreted in terms of an ellipsoidal model with axis ratios 
within several percent of a/b - 2.71 and b/c = 1.30 and a 
pole direction of either (71 ± 3°, +19 ± 3°) or (236 ± 
3°, +34 ± 3°), with the former direction strongly favored 
(Magnusson 1986, 1990). Drummond etal. (1991) obtained 
axis ratios of a/b = 2.56 ± 0.16 and b/c = 1.33 ± 0.01 and 
a pole direction within 9° of (69°, +10°), in reasonable 
agreement with Magnusson. A 250 x 90 X 70-km ellipsoid 
with Magnusson ’s favored pole direction would have pre- 
sented the same projected area to IRAS as a 135-km- 
diameter sphere; however, scaling such an elongated ellip- 
soid to the IRAS diameter in this manner is risky 
(Brown 1985). 

Occultation observations in 1980 and 1991 provide inde- 
pendent constraints on Kleopatra’s dimensions. The 1980 
occultation, near a time of lightcurve minimum (with S 
-15°), yielded five apparently complete chords from 40 to 
90 km in length plus four apparently incomplete chords 
that were fit by a 125 X 93 km ellipse (Dunham 1981). 
The 1991 occultation, near a time of lightcurve maximum 
(with S ^ -34°), yielded nine apparently complete chords 
fit by a 230 X 55-km ellipse (Dunham 1992), suggesting 
that the asteroid is even more elongated than lightcurve- 
based inferences would indicate. Nevertheless, for the mo- 
ment, let us adopt an ellipsoid with dimensions within 15% 
of 250 X 90 X 70 km and a pole direction within 3° of 
(71°, +19°) as an a priori model. 

The modest SNR of our Kleopatra echoes is adequate 
to reveal several noteworthy characteristics of the spectral 


signature. Figure 12 shows weighted sums of spectra in 
seven —45° phase intervals. The spectra from bins 4 and 
7 appear to be less than half as wide as those from the 
other bins, which is consistent with expectations about the 
asteroid’s extreme elongation. The shaded boxes define 
interval estimates of the spectral edge positions based on 
the a priori ellipsoid model. Since the phase origin of our 
spectra is not known accurately, we assume that the 
weighted average phase of the “group 4” spectrum coin- 
cides with the asteroid's minimum bandwidth orientation. 
Figure 13 shows two spectra obtained by taking weighted 
averages of the spectra from bins 4 and 7 (top) and from 
the other five bins (bottom). The shaded boxes in that 
figure correspond to the a priori model in its minimum 
and maximum bandwidth orientations. The boxes encom- 
pass the innermost zero crossings of the wide-bandwidth 
spectrum but are outside those of the narrow-bandwidth 
spectrum. To the degree that the zero crossings approxi- 
mate the true spectral edge positions, this result is consis- 
tent with the evidence from the 1991 occultation that Kleo- 
patra’s minimum pole-on breadth is shorter than that of 
the lightcurve-based model ellipsoid. 

The most intriguing aspect of the Kleopatra echoes is 
the bimodality of the wide-bandwidth spectra. As shown 
in Fig. 12 there is a central deficit of echo power in four 
of the five phase groups (1, 2, 3, and 5) that exclude the 
narrow-bandwidth orientations. The central deficit is seen 
most clearly in an average of all wide-bandwidth spectra 
(Fig. 13, bottom). This spectral signature, which resembles 
that of 4769 Castalia (Ostro et al 1990), is consistent with 
Weidenschilling’s (1980) conjecture that Kleopatra is a 
dumbbell-shaped asteroid. However, whereas high-SNR 
resolution of Castalia’s echoes in time delay as well as 
Doppler frequency and rotation phase allowed reliable 
estimation of that object’s shape and established its bifurca- 
tion (Hudson and Ostro 1994), our much sparser Kleopatra 
data set precludes shape reconstruction at a similar level 
of detail. Attempts to distinguish between a single and a 
binary shape for Kleopatra based on lightcurve data were 
inconclusive (Zappala et al 1983, Cellino et al 1985). 

Kleopatra’s weighted-mean OC radar cross section is 
7100 ± 1800 km 2 . The weighted-mean projected area of 
the a priori ellipsoid is 16200 ± 3400 km 2 , which is 13% 
larger than that of a sphere with Kleopatra’s radiometric 
diameter. (Here, weights based on the noise in each 4- 
min spectrum have been applied to the projected area at 
corresponding phases.) Our albedo estimate based on the 
a priori ellipsoid is & oc = 0.44 ± 0.15, where the assigned 
error is intended to include systematic uncertainties in both 
the radar cross section and the asteroid's dimensions and 
orientation. Kleopatra’s low “disc-integrated” circular po- 
larization ratio (p c = 0.00 ± 0.05) indicates that the radar 
echo must arise almost entirely from single back-reflections 
from surface units that are smooth at centimeter-to-meter 
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FIG. 12. Weighted sums of OC echo spectra of Kleopatra from 1985 within the rotation phase intervals depicted in the polar plot, each smoothed 
to a frequency resolution of 60 Hz. The frequency scale below spectrum 4 applies to all spectra. The shaded boxes show interval estimates of the 
spectral edge positions based on the a priori ellipsoid model {see text). 
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FIG. 13. Arecibo radar spectra of Kleopatra obtained in 1985. The panels represent averages over different ranges of asteroid rotation phase, 
as indicated by the insets (cf. Fig. 1). The solid and dotted lines plot the OC and SC echoes, respectively, which are smoothed to a frequency 
resolution of 60 Hz. The shaded boxes show interval estimates of the spectral edge positions based on the a priori ellipsoid model (see text). 
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scales. Under these circumstances, the OC radar albedo 
can be interpreted as the product gR of a gain factor and 
the Fresnel reflection coefficient at normal incidence (Eq. 
1). As noted in Section 3, the gain factor g would be unity 
for a smooth sphere and is expected to be within a few 
tens of percent of unity for most large main-belt asteroids; 
however, Kleopatra’s highly elongated shape could result 
in a gain factor that not only differs significantly from unity 
but also depends on the asteroid’s orientation. 

Discrepancies between Kleopatra’s radar signature and 
predictions of the a priori model suggest the potential for 
model refinement. Limited SNR confines this study to only 
the simplest models that can account for the spectra, but 
perhaps the most important benefit of this approach is 
that in modeling the echoes, one can estimate the average 
scattering law as well as the shape. In this way, the shape 
model accounts for structural contributions to g at “tile” 
scales, while the scattering law absorbs contributions to g 
from structure at smaller (“facet”) scales. This mitigates 
a significant source of uncertainty in the interpretation of 
Kleopatra’s radar albedo. 

A simple four-parameter model for Kleopatra consists 
of a biaxial ellipsoid (a ^ b = c) with a p * cos' 1 # scattering 
law. The model’s scale (in kilometers) is uncertain by at 
least —20% because of uncertainty in the pole direction; 
however, the scale factor cancels in the ratio a/b. As before, 
we assume that the weighted average phase of the “group 
4” spectrum in Fig. 12 coincides with the asteroid’s mini- 
mum bandwidth orientation. This has the effect of minimiz- 
ing estimates for the a/b ratio. Figure 14 shows the least- 
squares result, which yields a reduced x 2 of 1.10. The ex- 
pected correlation between the ellipsoid’s dimensions and 
the scattering law exponent was evident in our search for 
the best-fit values of a , b, p, and n; however, our estimate 
for this model’s a/b ratio, 4.1 ±1.3 (95% confidence), 
supports the hypothesis that Kleopatra is more elongated 
than Hghtcurve-based estimates. 

Two biaxial ellipsoids in contact provide a slightly more 
complicated model, with three additional free parameters: 
two for the second ellipsoid’s axes and one for the separa- 
tion of the ellipsoids’ centers. (The ellipsoids are allowed 
to overlap, but the a-axes are forced to be collinear. In 
addition, the Doppler frequency of the model’s center of 
mass is forced to coincide with the prediction ephemeris.) 
Figure 15 shows the least-squares result, which has a re- 
duced x 1 of 0.98. The —10% improvement in x 2 supports 
the hypothesis that the two-ellipsoid model provides a bet- 
ter fit to the data than does the one-ellipsoid model. The 
breadth ratio, defined to be the model’s total extent along 
the line containing the a-axes divided by the maximum 
extent orthogonal to that line, is 3.6 ± 1.1 (95% confi- 
dence). This supports our favored hypothesis that Kleopa- 
tra is more elongated than lightcurve-based estimates but 
does not rule out a breadth ratio as small as —2.5. 


The one- and two-component ellipsoid models provide 
acceptable fits to the data, but the 1991 occult ation chords 
(Fig. 16) suggest that both models are simplistic. Those 
chords, albeit few, show no evidence for a bifurcated shape. 
On the contrary, one anomalously long chord hints at the 
presence of a 25-km-high “mountain” (Dunham 1992). 
Extensive explorations of more complex shapes indicate 
that a variety of nonconvex shapes can provide good visual 
matches to our spectra but cannot survive statistical tests 
of significance — the SNR and geometric leverage of our 
data simply are inadequate. In light of the 1991 occultation, 
we cannot rule out the possibility that the bimodal spectra 
arise from radar reflectivity variations along the length of 
the asteroid. 

Regardless of whether or not Kleopatra is bifurcated, 
the ellipsoid models provide insight into the physical inter- 
pretation of Kleopatra’s radar albedo. It is useful to define 
the “equivalent spherical albedo” (<7 S ) to be the radar 
albedo that a sphere would have for any particular scatter- 
ing law (Hudson and Ostro 1994). The goal here is to 
separate the contribution of Kleopatra’s overall shape (at 
tile scales and larger) from the gain factor g, which may 
permit more useful comparisons with other radar-detected 
MBA’s, whose shapes are generally less exotic than Kleo- 
patra’s. For the cosine scattering law used above, a s = 
2 pl(n + 1). 

The equivalent spherical albedo of the one-ellipsoid 
model is fr s = 0-52 ± 0.21, and that of the two-ellipsoid 
model is &$ = 0.50 ± 0.19. Both estimates are similar to 
croc (the albedo based on the a priori ellipsoid), which 
shows that Kleopatra’s elongation need not result in a 
value of g significantly different from unity. On this basis, 
we expect physical inferences from comparison of Kleopa- 
tra’s albedo with other MBA values to be reliable. Further- 
more, if we assume that surface structures at facet scales 
can be treated statistically as a distribution of slopes, and 
if we adopt Parker’s (1973) slope probability density func- 
tion, then for a sphere with a p* cos' 1 # scatteri ng l aw we 
can write: g — 1 + SqI(2 ±^o), where s 0 = V2 Tn is the 
adirectional rms slope (see the Appendix). For n — 6, 
facet-scale surface structures increase a sphere’s gain factor 
g by only —14%. Thus, we expect Kleopatra’s radar albedo 
(froc) to be a reasonable first approximation to R. We will 
return to this issue in Section 5. 

12 Victoria 

Victoria’s lightcurves have a (synodic) period of 8.662 
hr (Erikson 1990) and a maximum amplitude of —0.3 mag 
(Lagerkvist el al 1989), similar to lightcurve amplitudes 
for Iris and Zelinda. Tempesti and Burchi (1969) estimated 
the pole direction to be (RA, Dec) = (4 h 12 m ± 6 m , +3 
±3°), noting that their quoted uncertainty is probably too 
small. The corresponding ecliptic coordinates are T1 = 



121 


radar 


OBSERVATIONS OF ASTEROIDS 



lot Phases are averaged 

the shaded boxes are as m F.g- 
ie (spectra 1-3)- 


l*r,*n; 4‘). A ‘^“'TmSel with M* 

,, 1W5) provides a ma ^ ^ , 1.00 ± 0.06 and ■ 
[OS of alb - t.25 coordinates of either D 

le direction with eel P _ 4 „ +40 £ 5 ). 

i 9*. +55 4=10-) “ “ ' ,Z e i al ellipsoid by noting 

^'D^d^po.. solutions yields 


and ft - -24 ' dLeter of 

“ t e We "'as a wording model for Victor,, 




122 


MITCHELL ET AL. 


(uncertain) 



occultation of 
19 January 1991 


FIG. 16 . Adapted from Dunham (1992). Chords from the 1991 Kleo- 
patra occultation (solid lines). Dashed lines are nondetections, one of 
which is uncertain. The longest chord suggests a large topographic feature. 


with a 15% uncertainty in each dimension, which is in- 
tended to include the quoted uncertainties in the radiomet- 
ric diameter and the ellipsoid's axis ratios, possible system- 
atic bias in scaling the ellipsoid using the radiometric 
diameter (Brown 1985), and uncertainty in the asteroid's 
orientation at the time of the IRAS observations. 

Victoria’s sky positions during our 1982 and 1989 obser- 
vations were ~30° apart. Figure 17 shows weighted sums 
of echo spectra for each apparition. Predicted ranges for 
Z? max based on the adopted a priori shape model and the 
Tl, Dl, and D2 pole directions are shown by the shaded 
boxes for comparison with the apparent spectral edges (the 
innermost zero crossings) for each year. As for Zelinda, 
we have used least-squares fits of an S(f) spectral model 
(Eq. 3) to place radar constraints on B. The open boxes 
in Fig. 17 show spectral edge positions from separate fits 
to each spectrum with n fixed to values between 1 and 6. 

Figure 18 shows constraints on Victoria’s pole direction 
based on our S(f) bandwidth estimates and the assumption 
that the maximum breadth sampled in each year was within 
15% of 132 km. For each year, our interval estimate for S 
defines an annular domain of possible pole directions (dot- 
ted curves). The intersection of the domains from 1982 
and 1989, outlined by the thick curves, bound the space 
of admissible pole directions. It is, of course, conceivable 
that an exotic shape and an unusually specular scattering 
law have conspired to produce spectra for which any of 


our bandwidth estimation approaches would lead to a gross 
underestimate of B\ however, to the extent that our experi- 
ence with the other targets discussed here is a valid guide, 
we prefer the hypothesis that our radar constraints favor 
the Dl pole direction over the D2 and Tl directions. 

Figures 19 and 20 show weighted sums of spectra within 
~50° phase intervals (a) and ^8° phase intervals (b) from 
the 1982 and 1989 observations, respectively. Figure 20 
shows that the 1989 spectra are persistently bimodal 
throughout the third phase interval, but are unimodal 
throughout the opposite (first) interval. This pattern is also 
evident in the somewhat noisier 1982 echoes (Fig. 19), but 
the rotation period is not known with sufficient accuracy 
to provide a common phase origin for both data sets. In- 
stead, we have numbered each of the 1989 phase groups 
to match the most nearly corresponding 1982 phase group 
based on the qualitative appearance of the spectra. 

We have explored the character of irregular shapes that 
can produce model spectra with Victoria’s unusual spectral 
signature; however, such models demand large numbers 
of shape parameters, precluding identification of a statisti- 
cally acceptable, demonstrably unique model. To the ex- 
tent that the backscattering is uniform and geometric 
(<x o (0) ~ cos 0), a Doppler spectrum maps the distribution 
of projected area across the object’s plane-of-sky projec- 
tion. An axisymmetric object spinning about a symmetry 
axis would give identical spectra at all phases. (Let us 
assume uniform scattering, homogeneous surface density, 
and principal-axis rotation throughout this discussion.) For 
a geometrically scattering target viewed equatorially, echo 
spectra 180° apart would be identical but with the Doppler 
frequency reversed. These extreme examples clearly do 
not apply to Victoria. Rather, the asteroid’s distinctive 
signature must involve some interesting shape and/or scat- 
tering properties. In light of (i) the persistence in phase 
of the spectral bifurcation, (ii) the apparently modest limb- 
darkening of MBAs at radar wavelengths, and (iii) the 
absence of observational evidence for (or theoretical argu- 
ments predicting) severe radar albedo variations on S-class 
MBAs, we expect that the object would have a strongly 
bimodal appearance in the orientation that produced our 
bimodal spectra. That is, we suggest that the simplest hy- 
pothesis is that Victoria’s strongly bimodal spectra most 
likely arise from a prominent concavity. Since opposite 
phases are not bimodal, we infer that the shape is not 
axisymmetric and that our view was at least a few tens of 
degrees from equatorial. 

5. RAMIFICATIONS FOR INDIVIDUAL OBJECTS AND 
ASTEROID SCIENCE 

Disk-Integrated Properties 

Table III lists for each experiment our estimates of the 
OC radar cross section, circular polarization ratio, and 
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FIG. 17. Arecibo radar spectra of Victoria obtained in 1982 (top) and 1989 (bottom) at the raw 10-Hz frequency resolution. Three lightcurve- 
based pole directions (Tl, Dl, and D2; see text) result in different predictions for the edge-to-edge bandwidth (stacked shaded boxes). The horizontal 
extents of the open boxes are determined by fitting each spectrum with an S(f ) model (Eq. 3) for 1 < n s 6. 
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FIG. 18. Comparison of radar and lightcurve constraints on Victoria’s pole direction shown in a rectangular projection of geocentric ecliptic 
coordinates. The radar constraints are obtained by fitting each spectrum shown in Fig. 17 with a spherical model that has a diameter within 15% 
of 132 km and a cosine scattering law exponent between 1 and 6. The aspect angle (a - <5 + 90°) between the radar line of sight and the asteroid’s 
spin vector is adjusted to maintain the best fit while the diameter and scattering law exponent are perturbed within the specified limits. This procedure 
yields an allowed range for a, which defines an annulus (dotted lines) centered on the asteroid’s coordinates (plus symbol). The sense of rotation 
is not determined with this method, so there are two annuli for each date corresponding to ± S. The intersections of the annuli (outlined in bold) 
define four possible regions for the pole coordinates that are consistent with both spectra. The radar constraints favor the Dl pole direction (A = 
9° ± 9°, /3 = +55° ± 10°) over the D2 and Tl directions (see text). 
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FIG. 19. (a) Weighted sums of OC echo spectra of Victoria obtained in 1982 within four —50° rotation phase intervals, smoothed to a frequency 

resolution of 30 Hz. Individual runs (or in some cases two runs that nearly co-incide in phase) are indicated by letters (a-c) within each group. 
The frequency scale below the group 1 spectrum applies to all spectra, (b) Weighted sums of OC echo spectra within each of the phase subgroups 
defined in (a). The weighted average phase (degrees) and phase subgroup are indicated, respectively, in the upper left and right of each spectrum. 
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FIG- 20. (a) Weighted sums of OC echo spectra of Victoria obtained in 1989 within three ~50° rotation phase intervals, smoothed to a frequency 

resolution of 30 Hz. Individual runs are indicated by letters (a-c) within each group. These are intended to correspond roughly to the group 
designations in Fig. 19. The frequency scale below the group 3 spectrum applies to all spectra, (b) The weighted sums of OC echo spectra within 
each of the phase subgroups defined in (a). The weighted average phase (degrees) and phase subgroup are indicated, respectively, in the upper left 
and right of each spectrum. 
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TABLE IV 

Average Radar Properties 


Target 

Class 

<Mc> 

{o’oc) 

Iris 

S 

0.14 (0.06) 

0.10 (0.03) 

Metis 

s 

0.14 (0.04) 

0.13 (0.03) 

Victoria 

s 

0.14 (0.03) 

0.22 (0.05) 

Kleopatra 

M 

0.00 (0.05) 

0.44 (0.15) 

Zelinda 

C 

0.13 (0.03) 

0.18 (0.06) 


Note. Weighted averages of fx c and & Q c from all Arecibo 
radar experiments for each target. For Iris, the uncertainty 
for may encompass real variations in fx c as a function 
of subradar latitude (see text). 


OC radar albedo, which were determined from weighted 
averages of all available spectra. For Iris, Metis, and Victo- 
ria, there is reasonable agreement in a oc from year to year 
once expected variations in projected area are accounted 
for. The circular polarization ratios for each target also 
show good agreement, except for Iris, as noted in Section 
4. Table IV lists for each target the weighted averages of 
/jl c and doc from all Arecibo observations. The low circular 
polarization ratios indicate that most of the echo power 
arises from single back-reflections from smooth surface 
elements, so that d- oc = gR may be used to constrain the 
normal reflectivity R. As discussed by OCS85 and Ostro 
et al (1991), for dry, particulate mixtures of meteoritic 
minerals with particle sizes no larger than —A/100, R de- 
pends almost linearly on bulk density (d), which in turn 
is a function of the pore fraction (porosity) and the zero- 
porosity density (the specific gravity) of the mixture's solid 
material. An empirical relation (Garvin et al 1985) that is 
applicable for porosities higher than —20% is 


d(R) = 3.2 In 


~ 1 + Vg ~ 
.1 - VR_ ’ 


( 4 ) 


For mixtures of metal and silicates with porosities lower 
than about 20%, small volume concentrations of metal can 
increase R by an amount that depends on the electrical 
properties of each phase and on the metal particles' dimen- 
sions and packing geometry. For solid enstatite chondrites, 
whose typical density corresponds to a reflectivity of 0.26 
by the above formula, laboratory investigations of “loaded 
dielectrics” predict reflectivities between 0.4 and 0.7 and 
measurements of meteorite specimens yield values of —0.7 
(Fig. 21b). 

Figure 21a plots the OC radar albedos of most radar- 
detected MBAs. When taken as a group, so as to average 
out target-to-target variations in g, S-class MBAs tend to 
have somewhat higher radar albedos than C-class MBAs. 
This is consistent with expectations about mineralogical 


differences between S- and C-class asteroids (Table V). 
The radar albedos of Iris and Metis are typical of radar- 
detected, S-class MBAs, which have an unweighted mean 
radar albedo of 0.14 ± 0.04 (OCS85). In contrast, Victoria 
has the highest radar albedo yet measured for an S-class 
MBA. Similarly, Zelinda's radar albedo is higher than most 
other C-class radar targets, which have an unweighted 
mean radar albedo of 0.11 ± 0.04. Differences in cr oc 
among our three S-class targets and between Zelinda and 
C-class targets as a group could be due to variations in 
surface mineralogy (e.g., Gaffey et al 1993), surface poros- 
ity, and regolith thickness; however, in the absence of de- 
tailed shape information, physical inferences based on the 
radar albedo of any individual target are obscured by un- 
certainty in g, especially in the case of Victoria, which 
evidently has a very unusual shape. 

Kleopatra's radar albedo exceeds those of all other ra- 
dar-detected MBAs (the next highest is 0.28 for 16 Psyche) 
and all radar-detected NEAs except 1986 DA (doc = 0.58; 
Ostro et al 1991). Moreover, the asteroid’s near-zero value 
of fi c distinguishes Kleopatra as one of the smoothest 
radar-detected targets at centimeter-to-meter scales. Can- 
didate meteorite analogs for M asteroids are irons (possibly 
with silicate inclusions) and enstatite chondrites, which are 
assemblages of NiFe metal and enstatite (Table V). With 
the above formula, a reflection coefficient equal to Kleopa- 
tra’s albedo corresponds to densities between 3.8 and 6.5 g 
cm -3 . Taken at face value, these numbers are consistent 
with a nearly entirely metallic composition and a porosity 
typical of the lunar regolith (—40%), but rule out an ensta- 
tite chondrite mineralogy unless the porosity is nearly zero. 
(Note that for any given porosity, R is two to three times 
higher for irons than for enstatite chondrites.) The equiva- 
lent spherical albedos derived in Section 4 argue that Kleo- 
patra’s high radar albedo is not simply an artifact of a 
highly elongated shape. Thus, unless the asteroid is unusual 
in some other respect (e.g., if the regolith thickness were 
m on Kleopatra but >1 m on other radar-detected as- 
teroids), it might be a core remnant of a collisionally dis- 
rupted, differentiated asteroid. 


Shapes 

Collisions are believed to be the dominant geologic pro- 
cess in the asteroid belt. Catastrophic impacts provide a 
natural explanation for asteroid families (Chapman et al 
1989), but serious questions remain about the extrapola- 
tion of laboratory impact experiments over many orders 
of magnitude to asteroid sizes and the extents to which 
internal material strength and self gravity play a role as a 
function of size (Fujiwara etal 1989 and references therein; 
Housen et al 1991; Ryan etal 1991; Nakamura et al 1992). 
Nevertheless, the shape distribution of fragments produced 
in laboratory impact experiments closely resembles that 
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FIG. 21. Asteroid radar aibedo versus d.ameter for selected MBAs. Each of our targets is identified by a 
by a symbol indicating taxonomic class. The radar albedcw of MB A tenets n^d^ii^wi m j^g^urvesThow R as"a function of porosity 

(b) Reproduced from Ostro er al (1991). Reflection coe cie / ^ 9g5 . Ulabv et aL 1990 ) over the range of lunar soil porosities, 
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meteorite specimens (Campbell and Ulrichs 1969). 


of moderate to small (D < 100 km) main-belt asteroids as 
inferred from their lightcurves (Capaccioni et al. 1984). 
Larger asteroids tend to have smaller lightcurve amplitudes 
and are thought to have more symmetrical shapes. Based 
on inferred collision rates in the main belt, large asteroids 
( D s 150 km) are thought to be thoroughly fractured by 
collisions (Farinella et al. 1982, Davis et al. 1989), possibly 
resulting in zero-strength -rubble piles’" with shapes gov- 
erned largely by the influences of gravity and centrifugal 
force (Weidenschilling 1980). However, little is known 
about the internal structures of MBAs, e.g., the size distri- 
bution of the inferred rubble and how internal structures 
might be manifested at the surface. An asteroid's ability 
to "reach an ideal equilibrium shape (e.g., a Maclaunn 
spheroid or a Jacobi ellipsoid) could be limited by the 
largest rubble fragments. 

Galileo images of Gaspra and Ida reveal concavities and 
flat regions spanning significant fractions of those asteroids 


dimensions (Thomas et al. 1994; Belton et al. 1994). One 
of the flat regions on Gaspra is 6 km across, or roughly one- 
third of the asteroid’s longest dimension. Stereographic 
measurements show that the maximum deviation from a 
plane in that region is 0.21 km. A preliminary analysis of 
Ida images reveals similar planar regions on an asteroid 
two times larger (Belton et al. 1994). The origin of these 
flat regions is currently unknown, although Greenberg 
et al. (1994) suggest that these regions and others could 
be of impact origin. In any event, such topography suggests 
the presence of coherent structures (possibly with signifi- 
cant tensile strength). Moore et al. (1994) suggest that 
stress-wave focusing could be responsible for grooved ter- 
rain on Ida, which would imply that the bulk of Ida is a 
single coherent structure. 

Radar spikes in our spectra of Iris, Metis, and Zelinda 
suggest that large flat regions exist on asteroids much larger 
than Gaspra and Ida— indeed, the flat regions themselves 
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TABLE V 

Asteroid -Meteorite Associations 


Meteorite properties 


Percentage 
abundance of 
Bulk metallics 

density 


Mineralogy 

Possible meteorite analogs 

(g cm -3 ) 

(wt) 

(vol) 

M asteroids 

Metal 

Irons 

7.6 

97 

97 

Metal + enstatite 

Enstatite chondrites 

3.6 

32 

19 

S asteroids 

Metal + olivine + pyroxene 

Stony irons 

4.9 

50 

30 

Metal + olivine + pyroxene 

Ordinary chondrites 
H 

3.6 

24 

12 


L 

3.5 

14 

8 


LL 

3.5 

2 

1 

C asteroids 

Hydrated silicates + 

Carbonaceous chondrites (CI1 and CM2) 

2.6 

11 

6 


carbon-organics-opaques 


Note. From Ostro et al. (1991). Meteorite analogs are from Lipschutz et al. (1989). Typical densities and approximate 
abundances of metallics (native metals and metal sulfides) are from Tables 2.1 and 4.1 of Dodd (1981). Table 24 of 
Buchwald (1975). Table 2 of Keil (1968). and Table 4.6 of Glass (1982). High-iron (H) and low-iron (L) chondrites are 
an order of magnitude more abundant than very low iron (LL) and enstatite (E) chondrites. The extremely rare Q asteroids 
might be analogous to ordinary chondrites (McFadden et al. 1989). 


appear to be larger than Gaspra. Are such features com- 
mon on asteroid surfaces? The fact that glints were ob- 
served from 3 of the 37 radar-detected MBAs may be 
significant considering the —10% probability of any partic- 
ular flat region being oriented within 10° of the radar- 
facing direction at some time during an observation with 
complete rotational coverage. Future observations (see be- 
low) are expected to provide a much larger and higher 
SNR asteroid radar sample. The prevalence of flat regions 
may improve our understanding of the internal structures 
of MBAs. 

Radar spectra of Kleopatra and Victoria suggest that 
the shapes of large MBAs can be far removed from simple 
ellipsoids, but the available information leaves us with 
more questions than answers. Kleopatra is by all accounts 
highly elongated, but it is unclear whether the asteroid is 
a “contact binary” or a single component with a prominent 
concavity or if strong radar reflectivity variations along the 
asteroid’s length have conspired to produce the bimodal 
radar spectra. Victoria’s unusual radar signature implies 
the presence of a prominent concavity, but the asteroid’s 
gross shape remains largely unknown. We note that the 
amplitudes and shapes of Victoria’s lightcurves are unre- 
markable (Lagerkvist et al 1987, 1988; Binzel et al 1989), 
which raises the question of what fraction of large MBAs 
possess large-scale topography and serves as a reminder 


of the limitations of triaxial ellipsoids as approximations 
to asteroid shapes. 

6. FUTURE OBSERVATIONS 

Our ability to constrain the shapes of Kleopatra and 
Victoria and the surface mineralogy of any of our targets 
was limited by the available SNR and the lack of delay- 
resolved echoes, which provide important geometric lever- 
age for the shape inversion technique described above. 
This situation is expected to improve dramatically upon 
completion of a major upgrade of the Arecibo telescope, 
which should provide a —20-fold increase in radar sensitiv- 
ity (Campbell et al 1994). Table VI lists OC SNR predic- 
tions for selected opportunities between 1997 and 2007 to 
observe our five targets with the upgraded Arecibo radar. 
(For comparison, the peak SNR per day was —2000 during 
the 1989 observations of Castalia.) These are by no means 
the only MBA opportunities — roughly five per year will 
have SNR > 100 per day. An experiment spanning several 
days on such a target would provide enough echo strength 
to place hundreds of delay-Doppler cells on the asteroid 
at each of many different rotation phases. Such data sets 
may be inverted to provide simultaneous estimates of the 
target’s three-dimensional shape and radar scattering prop- 
erties with the same techniques used for Castalia, thus 
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TABLE VI 

Future Radar Opportunities 


Target 

Current 
total SNR" 

Future 

SNR per day ft 

Year c 

Iris 

28 

100 

2000 



480 

2006 

Metis 

13 

140 

1997 

Victona 

54 

220 

2007 

Kleopatra 

19 

330 

1999 

Zelinda 

37 

1150 

2002 



1080 

2007 


"The OC SNR from the best experiment to date (see Table III). 
h Predicted OC SNR per day for nominal parameters of the upgraded 
Arecibo radar. For each prediction we assume an equatorial view and a 
radar cross section given by the product of the target's radar albedo 
(Table IV) and the projected area of a sphere with the target's radiometric 
diameter (except for Metis where we have assumed a diameter of 170 
km). We believe these predictions to be reliable within ~30%. 

c Selected opportunities between 1997 and 2007. The Arecibo upgrade 
is expected to be complete by mid* 1996. 


enabling us to place useful constraints on MBA mineralogy 
and potential meteoritic associations. 


Given the slope probability density function. p P (0), the adirectional rms 
slope (s Q ) can be calculated via (Simpson and Tyler 1982) 


^0 (tan $rms) 

= V2 in. 



tan : 0p P (0) sin OdB 


(A-5) 


Equations A-l through A-5 apply at the tile scale. A radar model is 
composed of a large number of such tiles, which collectively define the 
model's overall shape. As a concrete example, the radar cross section of 
a sphere with a diameter D and a cosine scattering law is 


~ I TT *dA = f (p cos" 0) sin Odd 
J dA 2 


‘ da 

JA' 

7 Tp 2 p 


2(h + 1) 

which corresponds to a radar albedo of 

- - 4<7 S 2 p _ n~ 2 

s nD 2 n + 1 n -r 1 

Finally, the sphere's gain factor g is 

frs n + 2 


S z 


R n + 1 


= 1 + 


2 + j?/ 


(A-6) 


(A-7) 


(A-8) 


APPENDIX 

Under the assumptions of geometric optics (facet > wavelength; 
Pc = 0) and perfect conductivity ( R = 1), Simpson and Tyler (1982) 
derive an expression for the radar cross section per unit surface area on 
the target, 

da sec 0 

= — — Pr(B)< (A-l) 

where p P (0) is the probability that a facet's normal is within a solid angle 
doj making an angle B w r ith respect to the mean surface normal (Parker 
1973). This probability density function is normalized via 

J" ~p P (0) sin Bd6 = 1. (A-2) 

To relate Eq. A-l to the empirical “cosine" scattering law (cr 0 = pcos n 0) 
introduced in Section 3 and utilized in our modeling of Kleopatra echoes, 
we first relax the assumption of perfect conductivity by multiplying the 
right hand side of Eq. A-l by the Fresnel reflection coefficient at normal 
incidence, R. Then, combining the cosine scattering law with Eqs. A-l 
and A-2 yields 

Pp(0) = (n + 2) cos n+i 6 (A-3) 


and 


R = f^. (A-4) 
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